A weighted network evolution model based on passenger behavior 
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This paper presents an evolution model of weighted networks in which the structural growth and 
weight dynamics are driven by human behavior, i.e. passenger route choice behavior. Transportation 
networks grow due to people's increasing travel demand and the pattern of growth is determined by 
their route choice behavior. In airline networks passengers often transfer from a third airport instead 
of flying directly to the destination, which contributes to the hubs formation and finally the scale- 
free statistical property. In this model we assume at each time step there emerges a new node with 
m travel destinations. Then the new node either connects destination directly with the probability 
p or transfers from a third node with the probability 1 — p. The analytical result shows degree 
and strength both obey power-law distribution with the exponent between 2.33 and 3 depending 
on p. The weights also obey power-law distribution. The clustering coefficient, degree assortatively 
coefficient and degree-strength correlation are all dependent on the probability p. This model can 
also be used in social networks. 

PACS numbers: 89.75.-k, 87.23. Ge, 89.40. Dd 



Most networks in real world are weighted. Recently 
weighted networks dynamics have attracted much atten- 
tion. People have developed many new models to un- 
derstand weighted networks evolution mechanism based 
on the Barabasi-Albert (BA) model [H which first pro- 
posed preferential attachment for networks evolution. 
Antal-Krapivsky (AK) model [2] relaxes degree prefer- 
ential attachment to weight driven attachment mecha- 
nism. Barrat-Barthelemy-Vespignani (BBV) model 0,0] 
adds local weight reinforcement dynamics to the struc- 
ture growth and weight-driven mechanism. Dorogovtsev- 
Mendes (DM) model gets similar result with BBV model. 
W. Wang et al[fij gives a global weight reinforcement 
mechanism to structure growth and weight driven attach- 
ment to model traffic networks evolution. Guimera and 
Amaral [2] models world-wide airport networks including 
spatial information as weight. Researchers begin to put 
more focus on weight dynamics and weight interaction 
with other factors [8j]. 

The above papers all assume the weight changes auto- 
matically at each step due to weight and topology cou- 
pling mechanism. However this assumption is too rough 
to understand underlying driving factors and helps lit- 
tle to control and manipulate networks. Transportation 
networks are good examples to explain this argument. 
Traffic is composed by passengers who are decision mak- 
ers. They decide which routes to choose and then changes 
the traffic on routes. The shortest path finding behav- 
ior and congestion phenomenon in networks have been 
studied by many papers 0, [l^, [HI, E2, Ell • Passengers 
collective route selection behavior play a critical role in 
the hubs formation and the pattern of network growth. 
For example, in airline industry airline companies pro- 
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vide low price for those who are willing to take transfer 
flights. Passengers make their decisions based on own 
time and money consideration. Some sensitive to price 
select transfer flights and contribute to traffic in hubs. If 
all passengers are not sensitive and don't want to trans- 
fer, the hubs and thus the popular hub-and-spoke oper- 
ation model in airline industry since 1970's cannot come 
into being. And the world-wide airline networks would 
not display the structure nowadays. So here we propose 
a new approach to understand transportation network 
evolution taking passenger's route selection behavior into 
consideration in addition to weight and topology coupling 
mechanism. 

The model assumes at each stage there emerges a new 
node with m destinations, i.e. new travel demand is cre- 
ated. A passenger who wants to go from origin to des- 
tination select his route. He can either fly directly to 
the destination or transfer at a third place. The proba- 
bility of taking non-stop flight is p and that of transfer 
flight is 1 — p. In the transfer case, the passenger deter- 
mines transfer node based on the weight between transfer 
node and destination. Analytical results show degree and 
strength both obey power-law distribution with the ex- 
ponent between 2.33 and 3 depending on p. The weights 
also obey power-law distribution. If p = 1, the model re- 
cover the strength-driven AK model. If we let 1 — p = S, 
it reproduces the results of the BBV model. Simulation 
confirms the theoretical results. This model is not limited 
in explaining transportation networks. It can also used 
to explain social networks which we will discuss later. 

The detailed model is described as follows. The initial 
network has No nodes with existing links. There are 
initial weights wq on the links. The networks evolutes 
according to the following rules: 

• Growth. In each step a new node n emerges with 
m destinations. The probability of one old node i 
to be the destination is determined by the ratio of 
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its strength Si to the total strength of all nodes: 



(1) 



which suggests that nodes with large strength have 
large possibility to be travel destinations. 

• Passenger's route choice. Then the new node n ei- 
ther connects node i with probability p or transfers 
at a third node j with probability 1 — p. In the for- 
mer case, a new edge with initial weight wq between 
n and i is established. In the latter case, things are 
little complicated. Passengers have to determine 
which node to choose as the transfer node. They 
only consider destination's neighbors which already 
have links with the destination. The probability of 
node jeu(i) being chosen as the transfer node is: 



Si 
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Here v(i) is the neighbors set of node i. Wij is the 
weight on edge between vertex i and j representing 
the traffic. Once node j is chosen, a new edge with 
initial weight wo between vertex n and j is created, 
and traffic between vertex i and j increases by wq. 
The quantity wq set the scale of weight. We can 
set wo = 1 without impact to the model. 

With the above two rules, the network grows and stops 
when nodes reach a preset large number TV. 

The rules has practical implications. In the first rule, 
a great number of traffic throughput indicate the city 
is more developed and more attractive to people than 
those city having little traffic. This rule is identical to 
those strength-driven models. In the second rule, there 
are two approaches to get to the destination. If passen- 
gers choose transfer mode, they prefer those nodes having 
frequent flights and large traffic between the destination 
and transfer node because frequent and large traffic can 
bring convenience, low price and less waiting time. 

This model is also applicable to social networks. For 
example, in author networks a new researcher is attracted 
to some famous scientist, but he has no approach to know 
him directly. Instead he finds that it's better to know 
the scientist's acquaintance with plenty of connections 
or cooperations between them. With the introduction 
of that acquaintance, he can finally get to know the sci- 
entist. This will also happen in friends networks, actor 
networks, etc. 

To analyze the distribution of degree and strength, 
we notice the model time is measured by the number 
of nodes added to the network, i.e., t = N — Nq. And the 
natural time scale of the model dynamics is the network 
size N. We treat strength S, degree K, time t as con- 
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tinuous variables 1 14] . When a new node n is added into 
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(a) transfer with X—p 



(b) connect directly with p 



FIG. 1: Visual description of evolution mechanism. New node 
n emerges at each step. The old node i becomes travel desti- 
nation with probability iSj/X^ Si. (a) Passengers select one 
neighbor j of node i, to transfer with probability Wij/Si. The 



weights update as Wij 



+ 1 and w n j = 1. (b) Passenger 



go to node i directly creating a new link between n and i. 
The weight Wni = 1. 



the network, the strength Si of an already present vertex 
i can increase by 1 if it is selected as a destination or 
increase by 2 if it is selected as a transfer node. And the 
degree K{ of the vertex i will only increase by 1 whether 
it is chosen as destination or transfer node. We have the 
following equations: 
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The equation J3]) can be rewritten if we notice the sum 
of strength at large time t is J2i Si ~ 2m(p + (1 — p)2)t 
and the initial condition Sj(i) = m,Ki(i) = m. So we 
obtain: 
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The strength and degree of vertices are thus related 
by a linear equation. It means the model is not only 
strength-driven but also degree-driven. 

From Eq. (jgj) we get the power-law distribution of 
the degree and strength P{K) ~ K~ r ,P(S) - S~ r with 
the exponent: r = 2 + 1/(3 - 2p) = (7 - 4p)/(3 - 2p). 
Obviously r is between 2.33 and 3 when p ranges from 
to 1. When p = 1, passenger always choose connect 
directly, or in other words there are no transfer behavior. 
The model is reduced to strength-driven AK model. If we 
set l—p = (5, the exponent becomes r = (3 + 4<5)/(l + 25) 
which reproduces the networks generated by the BBV 
model. It worths noticing because it indicates that a 



10° 10' 10 z 10 3 



10' 10 2 
K 



10' 



10' 
K 



10 4 



10 

w 



FIG. 2: Left: Probability distribution P(S). Data are consis- 
tent with a power-law behavior s~ r . The data are averaged 
by 20 networks of size N = 10000 and m = 2. In the inset, 
we give the value of r obtained by data fitting together with 
the analytical expression r = 2 + 1/(3 — 2p). Right: Correla- 
tions between strength and degree. The slopes are all around 
1 indicating linear relationship. 



model based on human behavior can lead to networks 
produced by automatic topology and weight coupling. S 
in BBV model represents the weight increase in every 
time period while 1—p represents the transfer behavior. 
The above finding indicates that transfer behavior can 
induce weight dynamics in BBV model. 

The time evolution of the weights Wij can also be com- 
puted analytically. Weight between already present 
node i and j evolves each time either j or i are selected 
as transfer node to another one. The evolution equation 
can be written as : 

dt m(1 p) \EiSi s t + E l s l sj (7j 

Noticing the initial condition Wij(Uj) — 1 and J2i &i ~ 
2m(p + (1 — p)2)t, we get 




(8) 



Thus the weights also have power-law distribution: 
P(w) ~ w~ a with the parameter a = (3 — 2p)/(l —p). 
When p = 0, a = 3. When the p — ► 1, a — ► oo. Actu- 
ally when p = 1 there are no transfer behavior and no 
weight dynamics. All weights are 1. So the probability 
distribution has an infinite parameter. 

In order to check the analytical predictions we per- 
formed numerical simulations of networks generated by 
using the present model with different probability p and 
network size N, The simulation recovers the analytical 



FIG. 3: Left: Probability distribution P(k). Data are consis- 
tent with a power-law behavior k~ r . In the inset, we give the 
value of r obtained by data fitting together with the analytical 
expression r = 2 + 1/(3 — 2p). Right: Probability distribution 
P(w). Data are consistent with a power-law behavior w~ a . 
We don't plot the curve for p = 1 because in that case all 
weights are 1 and the distribution cannot be displayed with 
an infinite exponent. The data are averaged by 20 networks 
of size N = 10000 and m = 2. 



results. The probability distribution of strength, degree 
and weight are in good agreement with theoretical pre- 
dictions. And the strength and degree have linear cor- 
relations. When we fit the data, we use the maximum 
likelihood method to estimate the parameter r See 
Fig © and Fig ©. 

Next we investigate the structural organization of the 
networks generated by our model by studying correla- 
tions of vertices. Using the formula defined by Eq. (4) 
and Eq. (6) of Ref [4-], we calculate degree-dependent 
cluster coefficient C{k) and average nearest-neighbor 
degree K nn {k) (also called degree-degree correlations). 
They both show disassortative behavior with C(k) and 
K nn {k) decreasing with k, which indicates there are hi- 
erarchal structures in the networks. The properties are 
depending on the parameter p. For large p, C(k) and 
K nn {k) are quite flat. While the p decreases, the curves 
grows. 

This disassortative behavior can be understood in the 
dynamical growth process. Vertices with large connec- 
tivities and strengths are the ones that enter the system 
early. New vertices are attracted to preexisting vertices 
with large strengths. In direct connection mode, there 
builds up a disassortative relation between "old" vertex 
with high strength and connectivities and "young" ver- 
tices with small connectivities. In transfer mode, the 
transfer vertices also have high connections because it 
generally is the most weighted in the neighbor of the 
destination. So there also builds up a disassortative rela- 
tion between transfer vertices and the new vertices. The 
average nearest-neighbor degree will be higher than that 
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FIG. 4: Top: Clustering coefficient C(k) depending on the pa- 
rameter p. Bottom: Average nearest-neighbor deg ree Knnik^ 
depending on the parameter p. The data are averaged by 20 
networks of size N = 10000 and m — 2. 



fer mode connects three vertices in one time period thus 
produces more correlations between vertices. With small 
p there are more transfer behavior in the networks lead- 
ing to large C{k) and K nn {k). 

In summary we have presented a model for weighted 
networks evolution that considers human behavior in ad- 
dition to weight dynamics and topology growth. We in- 
vestigated the evolution of degree, strength and weight. 
They are all distributed according to power laws with ex- 
ponents dependent on the probability p which determines 
the behavior of passengers. Clustering and correlations 
between vertices show clear disassortative behavior. The 
most difference between our model and previous models 
is that human behavior is incorporated to induce weight 
dynamics instead of automatic weight dynamics. This re- 
veals the underlying driving force in the growth of trans- 
portation networks and social networks. We believe that 
the model might provide a starting point for the realistic 
modeling that incorporates human behavior into techno- 
logical networks modeling. 



in direct connection because the transfer mode produces 
more highly connected vertices than the direct connec- 
tion mode. The cluster coefficient in transfer mode is 
larger than that in direct connection mode because trans- 
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